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The effective thermal diffusivity is evaluated for a two dimensional real solid and the real solid is 
modeled with periodic surface. The result contain the scale factor for the conformal transformation 
that flattens the surface. We find that the effective thermal diffusivity is always less than the local 
thermal diffusivity and this conclusion is independent of the model. 


Temperature and heat flow are two important quantities in the understanding the heat conduction. When temper¬ 
ature distribution is not uniform at all points of the body, then heat flows in the direction of decreasing temperature. 
The Fourier equation [Eq.(l)] is in general used in understanding heat conduction which is valid only for a stationary 
homogeneous isotropic solid without any heat source. But in reality we have solids of different shapes, e.g. we can 
have a cylinder with nonuniform cylindricity. Now the aim of this paper is to understand the effect of this type of 
nonuniformity on the thermal diffusivity. As an example we consider a two dimensional solid modeled by a peri¬ 
odic surface and we look for an analytical expression for thermal diffusivity. We start with Fourier equation of heat 
conduction [ij. 


where the constant a is called the thermal diffusivity of the medium, and is defined as 
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It combines three physical properties {k: thermal conductivity of the material, p: density of the material and Cp: 
specific heat) of the solid into a single constant. In Eq.(l) T{x,y,t) is the temperature at location {x,y) at time t. 
The differential equation of heat conduction will have numerous solutions unless a set of boundary conditions and an 
initial condition are fixed. We use standard no flux boundary condition 


n. V = 0 (on the boundary surface), 


( 3 ) 


where n is the unit vector normal to the surface boundary. The shape of the solid is defined by 

— oo < X < oo and yi{x) < y < y 2 {x). (4) 

yi{x) and y 2 (x) are the periodic functions of x with periodicity 27r, so that 

yi{x + 2 tt) = yi{x) i = l,2. (5) 


The effective thermal diffusivity can be evaluated by a similar procedure, as the one used in calculating the effec¬ 
tive mass of an electron under a periodic potential. In the following we use the same procedure, we look for the 
eigenfunctions of the following equation 

V' = -AV' (6) 

with n. v '0 = 0 at the boundary surface. As periodicity is involved in this problem, the eigenfunctions will have the 
popular Bloch form 


ifix+ 2TT,y) = e^^^^i}{x,y), (7) 

where k is the corresponding wave vector. Eor any value of k, there is a set of eigenvalues A„(fc). The smallest 
eigenvalue for any k is denoted by Ao(fc). So the effective thermal diffusivity a* is given by the limit 
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= a lim 
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Mathematically similar problems are already treated in a very interesting way by the use of complex variable technique 
ii using conformal mapping to flatten the surface boundaries 0] . We start with the conformal transformation from 
w = u + iv to z = X + iy. 


z = w + F{w), F{w+ 21^) = F{w). 


(9) 


By proper choice of -F, lines of constant v in the w planeare mapped into the surface boundaries in the z-plane. The 
period 0 < m < 27r maps onto 0 < cc < 2tt. The problem in using conformal mapping is the difficulty in finding the 
function F that flattens an arbitrary surface boundary. We follow the method used by Zwanzig [3] , chosing some trial 
functions F and comparing surface boundary it leads to with the one we are considering in this problem. Transforming 
from z-plane to the ic-plane we get the new Laplacian 


V ^ J 




where J is the Jacobian of the conformal mapping, and is given by 


J = 


dw 


dz 


So the eigenvalue equation in rc-plane becomes 


( 10 ) 


( 11 ) 


V = -jV” 

and ^ = 0 at n = ui and v = V 2 - The Block condition in w-plane will be 


The eigenvalue equation for (j) is given by 
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^ , d(j) , 2 , A 

I + 2ik— - k^(j) = 

ou J 


^ dv?- dv"^ _ 

Now integrating the above equation over u from 0 to 27r and over v from vi to V 2 , we get 


Cdv^du 


= \/k^ 
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Jvi Jo J(u,v) 
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It is natural to introduce 

a new periodic function (j), 



ijj{u,v) = 

(14) 

where 

4>(u + 27r, v) = 4>(u, v). 

(15) 


i:dvfrj^Hu,v) 

Now take the limit k ^ 0, the resulting expression of thermal diffusivity involves only geometrical scale factor, 


(16) 


(17) 


(18) 


In the following the shape of the solid is defined by its minimum and maximum widths i.e., ymin = 2[V — asinh{V)] 
and ymax = ‘2[V + asinh{V)]. For a solid with symmetrical shape, where V 2 = —vi = V, the effective thermal 
diffusivity is given by 


— = 1 + -a^ 


iinh{2V)\ 
2V J 


(19) 
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If we keep the minimum width constant and allow the maximum width to approach infinity, under that condition the 
effective diffusivity can be expressed as 


a* _ 4 

Umax 


( 20 ) 


The effective thermal diffusivity is always smaller than the local thermal diffusivity. For a solid of rectangular shape, 
the thermal diffusivity is independent of the width of the solid. But any deviation from the rectangular shape provide 
extra regions, that the ’heat’ can wander in before it can change its x-coordinate. So it takes more time to move 
in the x-direction, and hence thermal diffusivity is smaller. So our conclusion is true for any real solid and also for 
3-dimensional solid. 
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